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In this paper, we study the Loschmidt Echo (LE) of a two-dimensional Kitaev model residing 
on a honeycomb lattice which is chosen to be an environment that is coupled globally to a central 
spin. The decay of LE is highly influenced by the quantum criticality of the environmental spin 
model e.g., it shows a sharp dip close to the anisotropic quantum critical point (AQCP) of its phase 
diagram. The early time decay and the collapse and revival as a function of time at AQCP do also 
exhibit interesting scaling behavior with the system size which is verified numerically. It has also 
been observed that the LE stays vanishingly small throughout the gapless phase of the model. The 
above study has also been extended to the ID Kitaev model i.e. when one of the interaction terms 
vanishes. 



I. INTRODUCTION 



close to the AQCP and its collapse and revival with time. 



A quantum phase transition is a zero temperature tran- 
sition of a quantum many body system driven by a non- 
commuting term of the quantum Hamiltonian which is 
associated with a diverging length as well diverg- 
ing time scaleii^. In recent years, a plethora of studies 
are being carried out which attempt to bridge a con- 
nection between quantum phase transition and quantum 
information theorjs^i^. For example, information the- 
oritic measures like entanglement, quantum fidelity^^— , 
decoherence2r22 and quantum discor d^'^'^^ , etc., are be- 
ing studied close to the quantum critical point (QCP). 
These measures not only capture the singularities asso- 
ciated with the QCP but also show distinct scaling rela- 
tions which characterizes it. There have also been numer- 
ous studies on decoherence (or loss of phase information) 
in a quantum critical system which is closely connected 
to the LE to be discussed in this work; understanding 
decoherence is essential for successful achievement of the 
quantum computation. 

To study the LE in a quantum critical environment, 
we make resort to the central spin modeli^ in which a 
central spin S is coupled globally to an environmental 
spin model E (which in this case is the two dimensional 
Kitaev model). The LE (with the E in some ground state 
I ■00)) is given by 

L{t) = \{':Po\e'"°'e-'^"°+^"'^'\%)\^. 

Here the Hq and Hq + SHs are the two Hamiltonians with 
which the ground state |'0o) evolves, where the term SHg 
arises due to the coupling of the E with the S. It has 
been established that the LE shows a decay near the 
critical point of the E with a decay rate that marks the 
universality associated with the QCP of E^^—. Also the 
LE shows collapse and revival as a function of time when 
the E is at the QCP. 

The proposed work is organized in the following way: 
SecUpresents the model Hamiltonian, the phase diagram 
and discussion about the AQCP. In secHTl we describe 
the general calculation of the LE and in the subsequent 
subsections we study the scaling of the short time decay 



II. MODEL, PHASE DIAGRAM AND 
ANISOTROPIC QUANTUM CRITICAL POINT 
(AQCP) 

The Hamiltonian of the Kitaev model on a honeycomb 
lattice is given by 

j-\-l—even 

where j and / signify the column and row indices re- 
spectively of the honeycomb lattice while Ji, J2 and 
J3 are coupling parameters for the three bonds (see 
Fig. ([T])) ^^i^° ; and cr";, arc the Pauli spin matrices with 
a{= X, y and z), denoting the spin component. 
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FIG. 1: (Color online) Kitaev model on a honeycomb lattice 
with Ml and Af2 being spanning vectors of the lattice and Ji , 
J2 and J3, the coupling on the three bonds p^ . 

We will assume the parameters Ji, J2 and J3 are all 
positive and confine our analysis on the plane Ji + J2 + 
J3 = 4, since only the ratio of the coupling parameters 



2 



appear in the subsequent calcuations. The most excit- 
ing property of this model is that even in two dimen- 
sions it can be exactly solved using Jordan- Wigner (JW) 
transformationi2i2t2^ in terms of Majorana fermions 
given by 



6' 



Here, 
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(2) 



aj,i, a'- 1, bj^i and b'-i are all Majorana 
fermion operators, they obey the relations ^ = 

'^3,1^^1,1 = bj,i,{ajj,am,n} = {bjj,bm,n} = 2Sj^mSi,n and 
{o-j.iibm,i} ~ 0. One can now change the lattice site in- 
dices (j, I) of honeycomb lattice to a 2-dimensional vector 

n , where n = \/3ini + (^^i + f n2 which labels the 

midpoints of the vertical bonds of the honeycomb lattice. 
Here ni and n2 take all integer values so that the vectors 
n form a triangular lattice. The Majorana fermions 
and bfi are placed at the top and bottom sites respectively 
of the bond labeled by n. The whole lattice is spanned 

by the vectors Mi ~ — |j and M2 = + fj, see 
Fig. ©. 

Under the transformation to Majorana fermions as de- 
fined in Eqs. ([2]), Hamiltonian ([T]) takes the form 



H 



-Ml 



J2bnaf; 



JiDrAian) , (3) 



where D,-i=ib'fja'fj. These Dfi operators have eigenvalues 
±1 independently for each n and commute with each 
other and also with H which makes the Kitaev model ex- 
actly solvable. Since Dfi is a constant of motion one can 
use one of the eigenvalues ±1 for each ft in the Hamil- 
tonian. The ground state of the model corresponds to 
Dfi = 1 V n—. With Dfi = 1, we can easily diagonalize 
the Hamiltonian ([3]) quadratic in Majorana fermions. 

The Fourier transform of the Majorana fermions can 
be defined as 




e'*^-" + ale-'''-'' 

k 



(4) 



similarly 6^ also has same Fourier transform relation. 
The Ojt's and 6^'s are Dirac fermions which follow the 
fermionic anti-commutation relations. Here, N is the to- 
tal number of sites and N/2 is the number of unit cells. 
In the above sum given in Eq. (j?]), fc is extended over half 



of the Brillouin zone of the hexagonal lattice due to Ma- 
jorana nature of the fermions^S. We recall that the full 
Brillouin zone on the reciprocal lattice represents a rhom- 
bus with vertices {kx,ky) = (±27r-\/3, 0) and (0,±27r/3). 
In the momentum space the Hamiltonian ^ takes the 



ai,bi 

k' k 



and the re- 



form -ff = ^\^k'^k '^here ij^i = 
duced 2x2 Hamiltonian H-j^, can be expressed in terms 
of Pauli matrices as 



^k = "fc'^ 



-Pk<^^ 



where aj: — 2[Ji sin(fc.Afi) — J2 sin(fc.A/2)], 

and = 2[J3 -f Ji cos(fc.A?i) + J2 cos(A:.M2)]. (5) 

The eigenenergies of the are given by 



k 



■PI 



(6) 



This energy spectrum corresponds to two energy bands; 
it is noteworthy that for | Ji — J2I < J3 < (Ji -I- J2), the 
band gap A^: = — EZ vanishes for some particular k 
modes leading to the gapless phase of the Kitaev model. 
The phase diagram of the model is shown in an equi- 
lateral triangle satisfying the relation Ji + J2 + J3 = 4 
and Ji, J2, J3 > (see Fig. ([2])); one can easily show 
that the whole phase is divided into three gapped phases, 
separated by a gapless phase (inner equilateral triangle) 
which is bounded by gapless critical lines Ji = J2 + J3, 
J2 = J3 + Ji and J^ = Ji + J2. 

On the critical line J3 = J1 + J2 energy gap goes to zero 
for the four k modes given by {k^, ky)= (±27r/-\/3, 0) and 
(±27r/3,0) which are the four corner points of Brillouin 
zone. One can now expand and (3^ around one of the 
critical modes for J3 = Ji -I- j2, 



V3(J: 
3 



3(Ji 



Ji)kx 

{J, + J^)kl + -{.h + .h)kl + 



in the form 
J2)ky, 

2 , 3\/3 



{J2 



Ji)kxky, 
(7) 



where k^ and ky are the deviations from the above men- 
tioned critical modes. We note that varies linearly and 
Pj: varies quadratically in kx and ky. The point J^c = 2 Ji 
(where = J2) denoted by A in (Fig. ©) needs to be 
checked carefully. This is an AQCB^i with energy dis- 
persion Ej: ^ fc^ along fc^ {ky — 0) and Ej: ^ ky along ky 
{kx = 0). The corresponding dynamical exponents are 
given by zj^ = 1 and Z|| = 2, respectively 

For Ji ^ J2, Jsc = Ji + J2 is also an AQCP that can 
be shown using a rotation to a new coordinate system 



fci = \/3(J2 - Ji)kx + 3(Ji + J2)ky, 
k2 = 3(Ji -t- J2)kx - V3{J2 - Ji)ky, 
in which and fij: take the form 



(8) 



(9) 
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FIG. 2: (Color online) Figure shows phase diagram of the 
Kitaev model, satisfying Ji + J2 + J3 = 4. The inner equi- 
lateral triangle corresponds to the gapless phase in which 
the coupling parameters satisfies the relations Ji < J2 + J3, 
J2 < J3 + Ji and J3 < Ji + J2. Along the three paths I, II 
and III Ja is varied, so as to study the LE. The path I, II and 
III are defined by the equations Ji = J2, Ji = J2 + f and 
Ji + J3 = 4 respectively. 



'^JlJ2),C2 =27JiJ2(Jl + 

Ji) . Therefore, for a 



where ci = 9( Ji+ J2)(4 Jf +4 J2 
J2), and C3 = 18^/3JlJ2(J2 - 
general AQCP (Ji 7^ J2) the dispersion will vary lin- 
early and quadratically along fci and k2 directions, re- 
spectively, with two dynamical exponents zi — 1 and 
Z2 = 2. 



III. QUBIT COUPLED TO J3 TERM OF THE 
KITAEV HAMILTONIAN 



In this section we will provide a general calculation of 
the LE considering Kitaev model on a honeycomb lattice 
as an environment {E) that is coupled to a central spin- 
^ (S). We shall denote the ground state and excited 
state of the central spin S by \g) and |e) respectively. 5* 
is coupled to J3 term of E Hamiltonian only when the 
central spin is in the excited state |e). Therefore the 
composite Hamiltonian takes the form 



j-\-l—even 

+ J3(TliCTli+i+S\e){e\aliali^,), (10) 

where 5 is the coupling strength of 5 to iJ. We shall work 
in the limit of 5 — >■ 0. 

We consider that the S is initially in a generalized state 
l0(O))s = Cgl^) + Ce|e) (with the coefficients satisfying 
the condition |c„p + jcgp = 1), and the E is initially in 



the ground state |(p(J3,0)). The evolution of the envi- 
ronmental spin model splits into two branches, given by 
\ip{J3,t)) =exp(~iif(J3)t)|^(J3,0)) and |^(J3+<5,t)) = 
exp{—iH{J3 + S)t)\ip{J3,0)); the evolution of \'f{J3,t)) 
is driven by the Hamiltonian H{J3) — Ht{J3,0) (when 
the S is in the ground state and hence there is no S 
term present in the Hamiltonian), whereas |<p(J3 + S,t)) 
evolves with H{J3 + S) = Ht{J3,0) + Ve, where K = 
^Y^j+i=cvon'^j.i'^j.i+i' is the effective potential arising 
due to the coupling between S and E. The wave function 
of the composite system at a time t is given by 



m)) ^Cg\g)<E>\^{J3,t))+Ce\e)<S)\viJ3+S,t)). (11) 
As a result the LE is given by 

L(J3,i) = \{ipiJ3,t)\ipiJ3+S,t))f, 

= |(^(J3, 0)1 exp(-*iJ(J3 + 6)t)\^{J3,0))\\l2) 

Here, we have exploited the fact that the \(p{J3,0)) is an 
eigenstate of the Hamiltonian H{J3). 

Following Fourier transformation and Bogoliubov 
transformation the diagonalized form of the Hamiltonian 
dD) is given by 



HiJ3)^J2^-ej:{J3)AlAj: + ej:iJ3)BlBj:], 



(13) 



where the ^^r's and Sjt's are Bogoliubov fermionic oper- 
ators defined as 



1 



1 



V2' 



with 



"I 



^1 



"''hi 



(14) 



and the energy spectrum is given by (see Eqs.© and (jH])) 

^fe(J3) = + and e,r(J3 + S) = ^Jaj + Pf, (15) 

where aj: and /3j are defined in Eq.(l5]), and corre- 
sponds to the value with J3 + S instead of J3. 

The complete ground state of _ff ( J3) can be written in 
the form (sec Ref. [28| for details), 

i¥'(j3,o))=n[J(4 



k' ^ k k' 



(16) 



where k runs over half of the Brillouin zone of the hexag- 
onal lattice. Following mathematical steps identical to 
those in p"5lfl8| . it can be shown that Ea. ((T6|) leads to 
the expression for the LE given by 

i( J3 , t) = n ^fc = n [1 (20,^) sin^ (£,,( J3 +<5)t)] , (17) 



where, tan 6*^(73-^(5) = ot^l P'^ and 0^ = [Oj:{J3)-d^{J3 + 
5)]/2. The expression for LE closely resembles that of 
the case when the transverse Ising chain is chosen to be 
the cnvironmenti^. For numerical analysis of Eq. (|17p , we 
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shall use and ky in terms of two independent variables 



Vi and V2, with < Vi,V2 < 1. 



The fcj. and /cj, are given 



(18) 



which span the rhombus uniformly. Avoiding the corner 
points of the Brillouin zone ( where the LE results in a 
zero value) , we vary vi and V2 from 1 / {2N) to 1 — 1/ 
in steps oil/N, where N is the system size2^ and consider 
only the half of the Brillouin zone using the condition 

ivi,V2) > 1. 

The LE is calculated numerically as a function of J3 
using Eq.([T7|) and it shows dip at all critical points. To 
illustrate this, we choose three paths along which the 
interaction J3 is varied. In the first case J3 is varied 
along the path Ji = J2 (path T in Fig. ([2])) so that 
the model enters from the gapped phase to the gapless 
phase (extending in the region J3 G [0,2]) crossing the 
AQCP (point 'A' in Fig. @) at J3 = 2 - 5. The LE 
shows a sharp dip at point A and there is a revival with 
a small magnitude which again decays at the end point 
B, J3 = (see Fig. ([3])). Now surprising result shows up 
when the path is so chosen (path 'IF in Fig. given 
by the equation Ji = J2 + 1) that the system enters the 
gapless phase through an AQCP with Ji J2; (denoted 
by 'P' in Fig. The LE shows a sharp dip at the 

point P and stays close to its minimum value (with a 
small revival as observed in path I) throughout the gap- 
less phase and again shows a rise when the system exits 
the gapless phase through the point Q. In contrary, for 
the case when J3 is changed along the line Ji + J3 =4, 
J2 = (path 'III' in Fig. one observes only a sin- 
gle drop in the LE near J3 = 2 — (5 (see Fig. ([3]), inset 
(b)); this is associated with the critical point of the one- 
dimensional Kitaev model. In the next section we will 
study the scaling of the short time behaviour of LE close 
to these critical points and the collapse and reviaval of 
LE with time when the E is right at the critical point. 



A. Path I: Anisotropic Quantum Critical Point 

(Ji = J2) 

As discussed in SeclTTl J^c = 2Ji is an AQCP with 
critical exponents ly^ = zj_ = I along j direction and 

= 1/2, Z|| = 2 along i direction. At this point en- 
ergy gap vanishes for the three critical modes given by 
(27r/\/3,0) and (0,±27r/3) in half of the Brillouin zone. 
Now we will study the short time behavior of the LE (in 
Eq. (|17p ) close to the AQCP. We define a cutoff frequency 

= {K^ ci Ky such that modes up to this cut-off only 
are considered to calculate the decay of LE at short time 
close to the AQCP. Then the LE is given by 



Lc{h,t)= W 



(19) 




(a) 



(b): 
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FIG. 3: LE as a function of parameter J3 ( J3 is varied along 
path T) shows a sharp dip at point A ( J3 =2 — 5) and after 
a small revival in the gapless phase it again decays at ponit 
B, J3 = (see Fig. ^) with = Ny = 200, 5 = 0.01 
and t = 10. Inset (a) shows the variation in LE when the 
parameter J3 is varied along the path II (Ji = J2 4- 1) in 
the phase diagram for ~ Ny = 200, 5 = 0.01 and t = 10 
clearly showing a sharp dip at point P ( J3 — 2 — 5) and again 
rise at point Q (J3 — 1 — S). Inset (b) marks the dip in LE 
when J3 is varied along the path III (Ji -f J3 = 4), for this 
case A'' = 400, 5 = 0.01 and f = 10 so that E realise the 
change in the behavior at J3 = 2 — 5. Details of these three 
cases is provided in the subsection (|III All . (|III Bl) and (|III C|) 
respectively. 



We define the quantity S{J3,t), such that 5(J3,i) = 

In Lc = "'^k^ ky>o\'^^^k\- Expanding around 
one of the critical mode upto the cut-off, we get 
sin^ EkiJs + 5)t ^ 4{J3 + S- 2Jif and sin^ (20fe) sa 
/ {J3.-2Jif {.h + 5 -2Jif therefore we ob- 



tain. 



(J3 - 2Ji)' 



(20) 



where £ {Kc)^Att^Nc (iV^ + 1) {2Nc + 1) /5AN^ and is 
a integer nearest to 3NyKc/27T. We therefore find an 
exponential decay of the LE in the early time limit given 
by 



Le(J3,t) «exp (-7^2) 



(21) 



k^,k^>0 



where 7 = 36£ (K^) Jf6^/ ( J3 - 2Ji)^ The anisotropic 
nature of the quantum critical point is refiected in the 
fact that 7 scales as 1/Ny and is independent of Nx- 
Further using the expression of L^, {J^,t), one can eas- 
ily observe that it is invariant under the transformation 
Ny — >■ Nya, S ^ S/a and t — >■ ta, where a is some inte- 
ger. 

Now we fix Ji = J2 = 1 and J3 = 2 — 6 (point 'A' in 
Fig. ([2])) and observe collapse and revival of LE with time 
(presented in the Fig. The time period of collapse 

and revival is proportional to Ny, and is unaffected by 
the changes in N^ ; this confirms the scaling result of the 
decay rate 7 for the short time limit near the AQCP 
discussed above. 
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FIG. 5: (Color online) Variation of LE with t at the AQCP 
P (Ji — 3/2, J2 = 1/2 and J3 =2 — 5) shows collapse and 



revival with diffrent Nx 



Ny=Na.ndS = 0.001. 



FIG. 4: (Color online) The collapse and revival of LE with 
t at the AQCP (point 'A' in Fig. @) for Ji = J2 = 1, J = 
0.01 and J3 — 2 ~ S, keeping Nx{~ 100), fixed and varying 
Ny verifies the scaling relations satisfied by Ny, S and t as 
discussed in the text. The inset shows that the quasiperiod 
of the collapse and revival is independent of A^'^:. 



B. Path II: Anisotropic Quantum Critical Point 

(^1 7^ J2) 

It has been shown that the point P in Fig. ([2]) ( Ji 7^ J2, 
J^ c = J1 + J2) is an AQCP which can be seen by choosing 
directions ki ~ VS { J2 — Ji)i + S ( Ji + J2) j (see Sec |Tl| ^ 
and k2, perpendicular to fci^. The critical exponents as- ("^3,0 ~ exp (-7^ 

sociated with this critical point are given by i/i = zi = 1, 
and 1^2 = 1/2,^2 = 2 along ki and k2 directions, re- 
spectively. To calculate the early time scaling in a simi- 
lar spirit as in the previous section, we expand Eq. (|I7p 
near one of the critical modes up to the cut-ofF Kc to 
obtain sin^ (efe( J3 + S))t w 4 ( J3 + (5 - Ji - J2f and 



C. Path III: One-dimensional Quantum Critical 
Point (J2 = 0) 

As mentioned already, along the line Ji + = 4, 
(J2 = 0), the two dimensional spin model reduces to 
an equivalent one dimensional spin chain with energy 
gap vanishing at Ji = J3 for fee = tt and the corre- 
sponding dynamical exponent being z = 1. We shall 
now expand sin (efc( J3, 5)t) and sin (20^) around the crit- 
ical mode kc to analyze the short time decay of LE, re- 
sulting into sin^ (efe(^3 + 5)t) ~ 4 ( J3 + (5 — Jif P and 
sin^ (2(^fe) « Jfk^Sy (J3 - Jif iJ3 + 5- Jif. The LE 
hence takes the from 



(23) 



(20fc) « ki5^/A (J3 - Ji - J2r {j3 + 6^Ji^ J2r 



In the short time limit, the LE becomes 



Lc(J3,i)~exp(-7t2) 



(22) 




where, 7 = (S^^ (if,) / ( J3 - Ji - J2)', S {K^) = 
87r2j|iVc(A^c + l)(2A^c + l)/3iV^ and is an integer 
nearest to N Kd ^ttJ2- 

In fact comparing with the previous section fill Al one 
can see that in this case ki (instead of ky) appears in 
the expression of the LE in the short-time limit. Fur- 
ther, from Eq. and the expression of 7 one ob- 
serves that Lc{J3,t) is invariant under the transforma- 
tion Nx — Ny = N Na, 6 5/a and t — >■ ta, with a 
being some integer which is also observed in the collapse 
and revival behavior (see Fig. 



FIG. 6: (Color online) The LE as function of time at the 
QCP (point 'R' in Fig. (0) Ji = 2 and J3 = 2 - 5 with 
different N and 5 — 0.01, verifying the analytical scaling with 
N, S and t. 



where, 7 
.2 



4J^S^S{Kc)/{J3-Ji) and £ (Kc) = 
Att'^Nc (Nc + 1) {2Nc + 1) /6iV^ , Nc is an integer near- 
est to NKc/2tt. From Eq. we find that the LE 
shows a similar scaling relation as is expected for a one- 
dimensional chain with z = 1—; this is also confirmed by 
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studying the collapse and revival of LE (see Fig. ([6])). 

IV. CONCLUSION 

In this paper we study a variant of the central 
spin model in which a central spin (qubit) is globally 
coupled to an environment which is chosen to be a 
two-dimensional Kitaev model on a honeycomb lattice 
through the interaction term J3. Using the exact solv- 
ability of the Kitaev model, we have derived an exact 
expression of the LE when the interaction J3 is varied 
in a way such that the system enters the gapless phase 
crossing the AQCP of the phase diagram. However, the 
behavior of the LE as a function of J3 depends upon 
the path along which J3 is varied. In the case when the 
AQCP, Q (with Ji 7^ J2 see Fig. ([2])) is crossed, one ob- 
serves a complete revival of the echo when the system 
exits the gapless phase to re-enter the gapped phase; this 
is in contrast to the case Ji = J2. For the case of J2 = 
there is only one sharp dip at the critical point J3 = 2 — 5 
which is associated with the QCP of the one-dimensional 
Kitaev model. 

The early time scaling behavior for both the paths 
I and II close to the AQCP bear the signature of the 
fact that the gapless phase is entered crossing an AQCP 
with different exponents along different spatial direc- 
tions. This is also confirmed by studying the collapse 
and the revival of the LE as a function of time. However, 
one does not observe a perfect collapse and revival (ex- 
cept for the equivalent one dimensional case) ; this may be 
because of the proximity to a gapless phase. The quasi- 
period of collapse and revival in all cases scale with the 
system size as N^'. The case with J2 = reflects the 
fact that the system is essentially one-dimensional in this 
limit. It is straightforward to relate these results to the 
decoherence of the central spin close to a critical point. 

This study of LE can be verified experimentally as pre- 
sented by Zhang et ai^\ they measure the LE as an indi- 



cator of quantum criticality for a one-dimensional quan- 
tum Ising model with an antiferromagnetic interaction 
using NMR quantum simulators. In this experiment, 
they prepare the ground state of the Hamiltonian (using 
the gate sequences) which need not be the true ground 
state but could be a state that approximates the ground 
state of the system well, and then measure the LE for 
finite number of spins. Similar experiments can be re- 
alized with the approximate ground state of the Kitaev 
model. Also, since the Kitaev model can also be real- 
ized using an optical lattice^i^ (where the couplings can 
be separately tuned with the help of different microwave 
radiations), there exists a possibility of verifying these 
results in an optical lattice also. 

It should be noted that in a recent work, PoUmann 
et al^, have studied the problem of the LE in a trans- 
verse Ising spin chain in the presence of a longitudinal 
field; more precisely they calculated the magnitude of 
the overlap between the final state reached following a 
slow quench across the QCP and its time evolved coun- 
terpart at time t (generated following the time evolution 
with the final Hamiltonian). They observe a cusp-like 
minimum in the echo as a function of time in the limit 
when the spin chain is integrable. However, this behav- 
ior is smeared in the non-integrable case (with non-zero 
longitudinal field) thus providing a probe for integrable 
versus non-integrable behavior. In the present paper, we 
however deal with an equilibrium situation in which the 
spin chain is not quenched across the QCP, and observe 
the collapse and revival only at the QCP. 
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